Occupation numbers in Self Consistent RPA 

J. Dukelsky 

Instituto de Estructura de la Materia, Consejo Superior de 
Investigaciones Cientificas, 
Serrano 123, 28006 Madrid, Spain 
J. G. Hirsch 

0> ' Instituto de Ciencias Nucleares, Universidad Nacional 



X 



Autonoma de Mexico, 
Apdo. Postal 70-543 Mexico 04510 D.F. 



<: 

O ■ P- Schuck 

cn ■ 

Institut des Sciences Nucleaires, IN2P3-CNRS, 

' 53 Av. des Martyrs, 38026 Grenoble- Cedex, France 

00 ' 

■ A method is proposed which allows to calculate within the SCRPA theory 

o ■ 

On ■ the occupation numbers via the single particle Green function. This scheme 

complies with the Hugenholtz van Hove theorem. In an application to the 
^ ■ Lipkin model it is found that this prescription gives consistently better results 

C ■ than two other commonly used approximations: lowest order boson expansion 



and the number operator method. 
PACS numbers: 21.60.-n, 21.60.Jz, 21.60.Fw 

Keywords: Random Phase Approximation, Green function, occupation numbers, Lipkin 
model. 



I. INTRODUCTION 

The solution of the many body problem beyond the meanfield level is not a very well 
settled problem. Though the meanfield approach for all kinds of many body problems is quite 
uniquely defined, the determination of the higher order correlation functions is not. Besides 
the usual partial resummation of Feynman graphs ( e.g. ring summation in RPA) there also 
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exist variational ansatze such as those introduced by Jastrow, Gutzwiller, together with the 
Resonating Valence Bond approach, etc. [|I|. However only in rare cases these variational 
approaches can be worked through to the end by minimizing the groundstate energy so that 
any new route can have interesting perspectives. In most cases there remains the additional 
problem of how to determine the excited states. One of the attractive features of the Raleigh- 
Ritz variational Hartree-Fock (HF) theory is indeed that it yields, consistently within the 
same theory, groundstate and excited states (quasiparticle excitations). 

Since some time we have elaborated on a theory for two body correlations functions 
which in a certain sense can be considered as an extension of HF theory to two body clusters. 
We for instance obtain selfconsistent nonlinear equations for the correlation functions which 
simultaneously determine the correlated groundstate energy and the spectrum of excitations. 
We named this approach Self - Consistent Random Phase Approximation (SCRPA), since 
it is a consistent generalization of the standard linear RPA approach This formalism 

was also developed independently by a second group of authors which coined for it the name 
Cluster Hartree-Fock (CHF) which seems also very appropriate This type of theory took 
its roots several decades back starting with the work of Hara . Considerable progress was 
achieved by D. Rowe using the equation of motion method which is summarized in [H. Some 



years later the theory was rederived using the method of many body Green functions |110| , 
Since that time not much progress was made on the formal aspect of the theory until the 
more recent works cited above. 

The SCRPA has lately given a series of interesting results for various many-body prob- 
lems Nevertheless some open problems persisted in the past with this formalism 



concerning for instance the consistent evaluation of single particle quantities such as the 
single particle density matrix or the occupation numbers. An approximation which lately 
came very much in use in relating these quantities back to SCRPA (or to its poorer but 
numerically easier variant the so called Renormalized RPA (RRPA) [Q) is based on the par- 
ticle number method which long time ago already was advocated by D. J. Rowe W^. Very 



recently we have proposed and apphed a different method which calculates these quantities 



via the single particle Green's function with a mass operator coupling back to the SCRPA 
In those works, however, neither a detailed derivation nor an assessment of its quality 
was given . On the other hand it has been pointed out that certain consistency relations are 
indeed fulfilled. 

The purpose of the present paper is therefore to give a quite detailed derivation and to 
make a systematic investigation in a model case of the Green's function approach and to 
contrast it with other methods. 

The paper is organized as follows: In section II the SCRPA equations are deduced, their 
coupling with the single particle Green's functions is presented in section III, the application 
to the Lipkin model is developed in section IV, the numerical results in section V and the 
conclusions are given in section VI. 



II. OUTLINE OF THE PROBLEM 

Self Consistent RPA can be derived in various ways. The method which probably exhibits 
most clearly the analogy with ordinary HP theory is the one due to Baranger |T3|. Let us 
first rederive in this way single particle HP. To this end we define a mean single particle 
energy in the following way 

E I ^ - E^) I {N, 0| yp^a, \N + 1, u)f + W - E^') {N, 0| v^f^I W - 1, ^) " 



(1) 



E I (A^, 0| if^a, \N + 1, z/) r + {N, 0| v^fa^ \N - 1, u) 



where Ej; and |A^, u) are in principle exact eigenenergies and eigenstates of the Hamiltonian 
for a system with N particles. Por the groundstate we have u = and is a single particle 
creation operator. Minimizing (|I|) with respect to the amplitudes <y9^ and <y9^* leads directly 
to the following eigenvalue problem 



E 



H, al 



|0) = 



(2) 



where {...,...} is the anticommutator. 

It is easy to verify that (0) is just one of the forms of the usual single particle HF 
equations, once |0) is chosen to be a Slater determinant. 

Let us now in the same way find equations which describe another form of elementary 
excitations of the system such as density vibrations. To this purpose we define in analogy 
to (|I]) a mean excitation energy: 



u,k>k' 



E (i?.^ - E^) {N, 0| Xj:,,ala,, \N, v) - {E^^ - E^) {N, 0| Y,l,ala, \N, v) 



u,k>k 



{N,0\Xj:,,ala,,\N,u) - {N,0\Y,^,,alak\N, u) 



Minimization with respect to the amplitudes Xj^,,,Yt,, leads to 



(0| [SQ,[H,Ql]]\0) = E, (0| [SQ,Ql] |0) 



(3) 



(4) 



where 



QI — i-^kk''^k^k' '^kk''^k' 



k>k' 



a-k 



(5) 



and 5Q is a variation (with respect to X or Y) of Q'^ . Equation (^) constitutes the SCRPA 
equations which are described in great detail elsewhere Explicitly 
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(6) 



where the matrices A and B are double commutators coming from the left hand side of 
and M is the norm matrix to be discussed in the following section. They lead to a nonlinear 
eigenvalue problem for the amplitudes X and Y which therefore have to be determined 
iteratively very much like the HF eqs. (j^). Equations are equivalent to 



H,Ql 



\0) = E, 



(0| Ql H,Ql |0)=0 



(7) 
(8) 



This form is interesting since these equations have exactly the same structure as any 
mean field Hartree-Fock-Bogoliubov equations, be it for single Fermions or Bosons or, as 
here, for Fermion pairs. 

For a hamiltonian with two body interactions one verifies easily that (^) contains at most 
one and two body density matrices. Roughly speaking the two body density matrices can 
be expressed as quadratic forms of the amplitudes X and Y (for more details see 

An important point is to realize that (^) is not restricted to the particle-hole (ph) and 
(hp) subspaces as is of common use in the nuclear literature on the subject Here the 

only restriction in (^ is that it should not contain any diagonal (i.e. Hermitian) components. 
Therefore in Q\ = J2k=ik' Xkk'^^l^k' the matrix is not Hermitian. The single particle basis in 
which (§,||J^) shall be solved is obtained from 

(0|[i/,Qt]|o) = (0|[i7,g,]|0) = (9) 

One can show that (|^) is obtained from the minimization of the SCRPA ground state 
energy with respect to the basis but one also directly realizes that (^) is consistent 

with the equations of motion (1^,^. 

The matrix B contains the pair potential of the two fermion pairs whereas the matrix A 
contains the normal selfconsistent potential for Fermion pairs. Qualitatively we can represent 
the selfconsistent equations (||) as in Fig. 1 

Figure 1 

where the wiggly line stands for quantum fiuctuations. Such a selfconsistent mean field 
potential for density fiuctuations as shown in Fig. 1 seems quite natural, since the ground- 
state of an interacting Fermi system can be considered as a gas of quantal fiuctuations. The 
presence of fiuctuations also has a feedback on the single particle motion, an issue which we 
mainly want to consider in this paper. For example, to couple back consistently the single 
particle density matrix p^k' = (0| a^a^' |0) to the amplitudes X and Y in order to close the 
system of equations, has been a matter of debate in the past ||TB|. It should be noted that. 



depending on the problem at hand, it also can happen that certain elements of the two body 



density matrix can not directly be expressed via X and Y amplitudes. In the Lipkin model 
which we will study below we will see that indeed a particular matrix element of the two 
body density matrix falls into this category. We will, however, demonstrate that once we 
have a method at hand that allows to calculate the single particle density matrix we will 
also find a reliable method of how to evaluate the missing two body elements. 



III. COUPLING THE SINGLE PARTICLE GREEN'S FUNCTION TO THE SELF 

CONSISTENT RPA 

The eigenvalue problem (^) has as usual a corresponding Green's Function (GF) formu- 
lation. For the following it is useful to also briefly outline this approach which, of course, is 
completely equivalent to the eigenvalue problem (^. 

Let us therefore define the two time chronological Green's function at zero temperature 
which describe density fluctuations 

where T is the chronological operator and 

Oi = e'^'Oe-'"' (11) 

with H the full Hamiltonian operator. In principle in (|1^) one should take only the fluctu- 
ating operator a^a — (0|a^a|0) but since in the equations of motion (^,^ any c- number drops 
out we will stay with the definition given in (p!0|). 

The Dyson equation for ( ]T0| ) corresponding to reads after Fourier transformation in 
the approximation of the instantaneous mass operator [^^: 

kik2k^k2 kik2k-^k2 /-^ pxP2k^k2 ^ ' 

piP2kik2p-iP2 

with 

^k,k2k[k', = (0| [^l^kj^fk^ |0) (13) 



and 



(14) 



PlP2 



One easily recognizes from (|T2|-|I3p the equivalence with (^. Since the Eqs. (|T^,|TB|,|T^ 
have been derived at length in several preceding articles we will not represent them 
here. 

For the coupling with the single particle Green's function it is useful to define a SCRPA 
T-matrix from (0) in the following way 



G' 



r'O _i_ r'O rpsCRPA 

^kik2k[k'2 ^fcifc2PiP2 P1P2P2P2 P'lP2*^l'=2 



(15) 



with 



^0 

^kik2k[k!2 



nk2 - ^fci 



■^kik[^k2k'2 



(16) 



where = (OlalafclO) and Ek 



iO - El + 62 

^kk'kk' '^fc' occupation numbers and 

generalized single particle energies which we assumed without loss of generality to be diag- 
onal and Vkik2k3k4, is the antisymmetrised matrix element of the two body interaction. With 



r^-[T^) the T-matrix in ([15|) is uniquely defined. Since this is quite standard procedure we 



do not further elaborate on the form of the T-matrix. A form equivalent to (|T5|) is given by 
(we use summation convention) 



r< _ r^o I ^0 T^SCRPA 

<-^fcifc2fcifc^ - '^kik2k[k'2 + <-^fcifc2PlP2^PlP2P^,P^ *-^PiP2'=ifc2 



(17) 



with 



T^SCRPA 
^kik2k[k'2 



^k!k2k[k'2 " i^ki - ek2) hik'^hik'^ 



(18) 



From (|T7|-0) we also read off the equality 



^kik2kik4,^k-ikAk[k'2 — _^-'fclfc2fc3fe4 '~^k3k4k[k'2 



(19) 



^nSCRPA ryO 

The important point to recognize is that the mass operator of the single particle Dyson 
equation 



iu^-e,) G-^,=6,,,+J2Mt,G;,, (20) 

p 

has a well known representation in terms of the full two body T-matrix For better 
visibility we present the relation graphically in figure 2. 
Figure 2 

At this point it has now become obvious what our interrelation of single particle GF 
and SCRPA shall be: we have to replace in Fig. 2 the full T-matrix by the approximate 
rpSCRPA defined in ([T3|). In addition to this obvious construct there also exists a direct 
and strong consistency requirement. It stems from the fact that we have now two ways of 
calculating the correlation energy: the first uses the well known relation between the single 
particle GF and the ground state energy |^,|16|,|17 



Eo = -'- lim Trf^4 + .,) Gl/ (21) 



2 t'-t^o+ \ dt 

The second expresses the correlation energy density via the two body GF (|T0|) : 



Ecorr = 7 lim Tr 
4 t'-t^o+ 



where again Vkik'^k2k[ is the antisymmetrised two-body matrix element entering in the Hamil- 
tonian H . 

The requirement is now that both expressions for the correlation energy, that is, the one 
deduced from (^) and (^2]), agree. This is equivalent to the Hugenholtz-van Hove theorem 
which states that the chemical potential calculated via the single particle GF must be 
equal (at equilibrium ) to the energy per particle when calculated from (^21). It turns out 
that this only is achieved when expanding the GF in (pO|) to first order in the mass operator 

G, = Gl + Gl Mt Gl (23) 

with 

{u - e,) Gl = 1 (24) 
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Of course one can use the iterated solution of the Dyson equation, i.e. Gk = 
{u — Sk — but for consistency then the particle-hole propagators of the SCRPA equa- 

tion must be redefined accordingly. This has been discussed in and may be elaborated in 
the future but for the moment we keep with the more restrictive consistency relation ( |23D 
together with (|T^-|T5|). 

For space reasons we have been relatively short in this general section. We will, however, 
work out in some detail the model case of the next section so that the reader, by analogy, 
shall be able to reconstruct details also in the general case quite easily. 



IV. APPLICATION TO THE LIPKIN MODEL 



The Hamiltonian of the Lipkin [0 model is given by 

V 



H = eJo-^{Jl + J^_) (25) 



with 



m=l 

n a 

J+ = ^ c\j^Com , J- = ^ cl^Cim (26) 

m=l m=l 

The indices and 1 denote the lower and upper levels respectively, separated by an energy 
e, and m is the angular momentum projection in each shell with degeneracy Q. 

The commutation relations between these three operators, which are the generators of 
the SU(2) group, are 

[J+,J_]=2Jo, [Jo,J±] = ±J± . (27) 

In the Lipkin model the number of particles is exactly that needed to completely fill the 
lower shell, i.e. N = Q. 
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A. SCRPA equations 



The SCRPA solutions are built with the operators (we stay in the normal phase) 



1 



1 



acting over a correlated vacuum |0), which is defined by the equation 

Q\RPA) = 



(28) 



(29) 



to yield the excited state 



\1) = Q^RPA) 



(30) 



The SCRPA equations (^ then take the following form 



A B 
-B -A 



X 
Y 



E 



X 
Y 



with the matrix elements A and B defined by fT^Jig 



A={[J_,[H,.U^]])I{[J_,.U]) 
i?=([J+,[//,J+]])/([J_,J+]) 

where we used (■ ■ ■) for {RPA\ ■ ■ ■ \RPA) . 

The normalization of the excited state Q'^\RPA) is given by 

m^) = {[Q.Q^]) = x^-y' = ^ 



(31) 



(32) 



(33) 



With the inversion of (|2S|) yields J+ = ^-2{Jq) [XQ^ + YQ) and the matrix 
elements of the SCRPA matrices read 



AscRPA = e + 2VXY 



B 



SCRPA = 2V^-/fl + V{X' + Y') 
wo/ 



(34) 
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From (0) we see that we face exactly the problem discussed in Sect. 3. The single 
particle occupation (Jo) and the square (Jq) can not directly be expressed in terms of X 
and Y . Of course as is well known in the present simple model it is possible to calculate 
the RPA groundstate via (^) explicitly ■ 

I«^^>-|^|^^(^)^4'I^^) (35) 

and therefore also (Jo) and (Jq) can explicitly be calculated However, this is not the 
usual situation and in general it will be very difficult if not impossible to solve the vacuum 
condition We have therefore to develop other methods to get access to these quantities 
independently. As a word of caution we should mention again (see Section II) that it is 
not possible to include Jq as a further component into the definition of the RPA excitation 
operator Q'^ (|28|) since Jq is hermitian and it is then impossible to define the norm of the 
RPA excited state. 

In the next section we will therefore elaborate on the evaluation of (Jq) via the single 
particle GF the way we have outlined it in the general section 3. 

For later use we also introduce here the renormalized RPA (RRPA) matrix elements 

AnRPA = e, BRnpA = 2V{Jo) . (36) 

They will be used below when we will compare the results of the SCRPA not only to the 
exact solution but also to RRPA. 



B. SCRPA and the single particle Green's function 



As outlined above we have to construct a mass operator for the s. p. GF such that it 
yields exactly the same groundstate energy via eq. ( PT| ) as when calculated directly from the 
two body GF (^). In order to explain the principle we first want to exemplify the procedure 
with standard RPA. In this case we have to put in eq. (0) X = 1 , Y = and (Jq) = — ^, 
(^o) ~ X ■ ^'^^ example consider the interaction energy to RPA order 
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E,ot = ~ ((4) + ( E^^^ = -VnXY (37) 
Using one of the RPA equations (dropping 1/Q corrections): 

V nX = {E + 6) Y (38) 
and multiplying this equation with X we obtain for Epot ■ 



E^^t^ = -Vn^jn (39) 



Expression ( |55D can be identified with the evaluation of the Feynman graph shown in 
Figure 3 
Figure 3 

where the wiggly line represents the RPA phonon with energy E . The particle-hole bubble 
has energy Sph = e and together with the phonon the vertical cut has energy E + e what 
corresponds to the energy denominator in (|39D . The amplitude of the phonon is and the 
two dots of the graph represents the interaction squared. As usual we can obtain the mass 
operator from the groundstate graph in cutting open the hole line. Therefore we obtain e. 
g. for the GF of the upper level ( G^^* = —i (t (^aim (t) a[,^ (^'))) ) approximation 
ofeq. (0) 

OJ— 2 ^ ~ 2 00 + rj + ^ ~ 2 

where the mass operator has the obvious graphical representation of Fig. 4 
Figure 4 

Using the (exact) relation (what is just a variant of (0)): 

1 lim (^^-^]G\-;^' = -}^/jl) (41) 

2 t'-t-.o+ V <9t 27 ^" 2 \ +/ ^ ^ 

and inserting into the Ihs expression (^) we obtain 

(^l-^]G\-^' = J-(vnf^ (42) 
t'-t^o+ \ dt 2J 2^ ' E + 6 ^ ^ 
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This is just half the potential energy (^) in the standard RPA approach. Proceeding 
analogously with Gom and adding to ( |^ ) the corresponding expression yields the missing 
factor 2 . This demonstrates that our construction of the mass operator in ( ^OD is consistent 
with the RPA groundstate energy. 

Since we now have the s. p. GF at hand it is straightforward to calculate the occupation 
numbers via 

aL^im ) = , lim G*!"*' (43) 



Inserting into (^31) the rhs of (^OD yields 

E = (^^)' T^f^ = (44) 

m \ ' J 

where in the last equality we again made use of the RPA equations. It is now easy to restore 
the value for (Jq) since ('^om'^om) = ^- Y.m ('^L'^im / and therefore 



(^o) = -^ + y' (45) 
It is interesting to realize that (^) corresponds to the Holstein Primakoff boson expansion 



of (Jo) [|T4|; a result which of course is consistent with RPA theory. 

Let us now repeat the same procedure but with SCRPA. Using the SCRPA equations, 
in analogy to the steps above, we can write for Epot '■ 



El^r' = -2 (Jo) V^^^^^f^ (46) 



where A = A — e and A, B are determined in (0). Again in cutting open the hole line we 
now find in analogy with ( ^OD for the mass operator according to (p3|) 



'"^ oj-l cu-| uj + E + l 



(0|J_|1)A+(1| J+|0)i? (47) 



where |1) again is the excited state |0) . In the RPA limit we obtain (^Ol). We immediately 
check that indeed we get back from Gim ( and Gom ) the correct expression ( ^BD for E^^^^^ 
inserting ( ^T]) into the Ihs of (|l|) ( and similar for Gom )• 
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Since we now have a consistent SCRPA expression for the single particle GF at hand we 
proceed, as this was our goal, to the calculation of (Jq)- Inserting ( ^Tj) into the rhs of (^31) 
one directly obtains 

n 
n 

2 



1+2XY-, 



(e+E) 

Of course this is still an implicit equation for ( Jq), since the SCRPA eigenvalue E depends 
on it. Before proceeding it is interesting to study several limits of (48). Of course for the 



interaction going to zero we recover the free gas limit (Jq) = ~^ ■ We already checked that 
(^) goes over into the RPA limit ( ^0|) when A, B and the transition amplitudes are replaced 
by their RPA expressions. Therefore we also recover the boson expansion result. 

One should note that in order to obtain the correct RPA result one must not make the 
mistake to go over to the RPA limit, i. e. X = 1 , F = , (Jg) = (Jo)^ = x 
(^) because to get (^8]) already the assumption has been used that (Jq) 7^ — f on the rhs 
of (^) what would not be consistent with the RPA groundstate energy then. Now if we 
nevertheless take the RPA limit, using directly one obtains 

W - -y^yi • («) 

This result is interesting because it is precisely the lowest order result which one obtains 
with the number operator method []T3|,^ . In the light of our theory this formula ( P^D seems 



to be inconsistent because if on the rhs of (^) one keeps (Jo) 7^ — f , there is no reason to 
drop all the other terms going beyond standard RPA. So in this light the pure lowest order 
boson result ( ^5]) seems to be more consistent than the partially resummed series (^91) . We 
will see later that this is indeed confirmed by numerical results. 



C. Determination of (Jq J 



In principle we are still short of the expectation value of the square of the occupation 
number. Eventually we could try to establish an analogous expression to what has been 
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found for (Jq) (0). However, at least in the present model the factorization relation 

(^'> = (^o)' (50) 

seems to be extremely well fulfilled for the whole range of the interaction strength considered 
(see next section). Of course this may be a particularity of the model but we suppose that, as 
long as the operator Jq or analogous operators in other problems are sufficiently collective, 
equation (^) should work quite reasonably. In order to check this we present the ratio 
^ = —a/ {Jq)/ (•^) for the fixed interaction strength x = y (fi — 1) /e = 1. (i. e. at the 
meanfield transition point where fiuctuations are expected to be maximal) as a function of 
Q in Figure 5. The exact results are represented by full squares, those obtained using the 
exact RPA vacuum (|5D by a full line (SCRPA) and a dotted line (RRPA). 
Figure 5 

Only for Q values lower than 4 one can see a significant deviation from unity. So definitely 
s-wave shells are difficult candidates. On the other hand should there be no degeneracy at 
all like in a rotating nucleus or in an electron system in a magnetic field there is no need to 
know the occupation number square since we have anyway 

ai'^koi'^k) = \a\^k) (51) 



So unless there is appearance of two fold degenerate levels in a problem one is prob- 
ably well off with the factorization (|50|) . In the former case a perturbative expansion of 
square operators in terms of linear operators as proposed in [|| using RPA excited states 
as intermediate states should adequately improve on (|50|) which represents the zero order 
approximation. This approximation is based on expanding the expectation value of any two 
body operator by inserting a complete set of RPA states. Specifically for the Lipkin model 
we have 

Truncating to first order and evaluating the expectation values using the vacuum condi- 
tion we finally arrive to 
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a relation which expresses ( Jg) in terms of (Jo) • 

Let us next study the numerical results as they follow from our SCRPA theory described 
above. 

V. NUMERICAL RESULTS 

In this section we mostly will present results for Q = 14. We will begin in first place to 
investigate the quality of the results for the correlation part of the groundstate energy, i. e. 
the correlation energy 

E^rr = {H) - e- (54) 

with 

(H) = (Jo) [e-VXY] (55) 

Figure 6 

We show Ecorr as a function of x = ^ ^ 1) /^^ in Fig- 6 for the RPA (dashed line) , 
RRPA (small dots), SCRPA (full line) and the exact solution (full squares). It is a very well 
known fact that RPA due to the quasiboson approximation, i. e. the violation of the Pauli 
principle, overestimates in general quite strongly the correlations and in fact overbinds in 
the groundstate energy. This the more, the closer one comes to the phase transition point 
where RPA collapses. This strong overbinding of the RPA was for example also found in a 
recent calculation of the electronic binding energy of a metallic cluster. When compared 
with the exact results SCRPA performs extremely well for E^orr up to and even beyond the 
mean field phase transition point x = 1 whereas RRPA starts to deviate strongly from the 
exact result at x = 1 . 

Since it is not possible to distinguish in Fig. 6 that the SCRPA values of Ecorr stays 
consistently above the exact ones, we also present the results in Table 1. 
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Table 1 

Another interesting quantity is the excitation energy. We show ii^ as a function of x 
Fig. 7. A similar scenario as in the previous figure prevails: SCRPA yields by far the best 
agreement with the exact results though the differences for x ^ 1 ^^e now more pronounced. 
It is also true here that the SCRPA excitation energy stays consistently above the exact 
results as can be seen from Table 2. 

Figure 7 

One could conclude from that that the SCRPA also leads to an upper bound for the 
excitation energy. This conjecture may be backed from the fact that we actually derived in 
Section 2 the SCRPA equations from a minimization within respect an average excitation 
energy. However, before drawing any definite conclusion in this respect, a more general 
model with more levels must be studied. 

Table 2 

Let us now come to the investigation of the quality of the different expressions for (Jo)- 
There are essentially three: the one which we prefer on theoretical grounds is the one from 
the Green's function approach since it is the only one which fulfills a strong consistency 
relation with SCRPA equations (i. e. the Hugenhotlz-van Hove theorem). The second is the 
quasiboson approximation (^) which represents the lowest order correction in to the 
free gas results. The third comes from the so-called number operator expression (^) which 
has recently become very popular in the nuclear physics literature |[T5|j20|j2^ . We have shown 



that it is as well obtainable from the GF approach in operating additional approximations 
to (^HD, and that those approximations are not consistent among them. 
Figure 8 

In Fig. 8 we show the quantity VL/2 + (Jq) as a function of x for the three approximations 
to (Jo) when used in the SCRPA equations (of course only the one corresponding to GF 
method corresponds to our definition proper of SCRPA). In addition we show in Fig. 8 also 
the exact result (full squares). The solution of the GF method (and therefore the SCRPA 
proper) is shown by the full line. The quasiboson approximation is shown by the broken line 
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and the number operator method by the dotted hne. Not unexpectedly the GF results are 
closest to the exact ones. Somewhat a surprise is that the number operator method works no 
better than the quasiboson approximation. However in the light of our discussion in section 
4 where we argue that one passes from the GF expression (^) in an essentially uncontrolled 
way to the number operator expression (^Sp this outcome may seem less astonishing. 

We should also say that the injection of (Jo) and ( Jq) as expressed with the RPA ground- 
state wavefunction (^) into the SCRPA equation still improves the results in Fig. 8 with 
respect to GF. However, we do not show this result in order not to overload the figure and 
because it corresponds to a situation which in general is not realizable. 

In Figure 9 VL/2 + (Jq) is shown not as a function of x for fixed Q but for fixed x = 1- 
function of f2 . 

Figure 9 

Again we see that Q = 2 appears as the worst case. It is, however, interesting to see that 
for this case the differences between the various approximations are also largely enhanced 
without, however, inverting their respective order. 

One last interesting quantity is the ratio Y/X as a function of x, shown in Fig. 10. It is 
well known that this ratio goes to 1 when approaching the phase transition point in RPA (as 
seen in th broken line) while the value of X and Y tend to oo individually. This then makes 
any RPA result close to a phase transition meaningless. On the other hand in SCRPA this 
ratio still stays of the order 1/2 around the transition point and also X and Y remain within 
very reasonable limits (X = 1.156, Y = 0.580 at x = 1.) 

Figure 10 

A word of caution is worth here. While the energetics and the occupation numbers 
obtained with the SCRPA are very close to the exact ones, the wave functions around and 
beyond x = 1 (the value at which standard RPA collapses), being far better than those 
obtained with RPA or RRPA, can nontheless have an overlap with the exact wave function 
of less than 50% P3. In this case the SCRPA must be extended to the deformed basis IBIl. 
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VI. CONCLUSIONS 



In this work we addressed the question of how to close the SCRPA equations in a consis- 
tent way and, in particular, of how to calculate single particle quantities such as occupation 
numbers in this formalism. We showed in detail how to couple back SCRPA into the single 
particle propagator consistently. The consistency criterion was based on the fulfillment of 
the Hugenholtz-van Hove theorem which states that the chemical potential obtained from 
the single particle propagator must be equal (at equilibrium) to the energy per particle 
when directly calculated via the correlation function. For some problems (for instance in 
such schematic models as considered here) there may also be correlation functions which 
involve the expectation value of the square of the occupation number operator, which fall 
out of the SCRPA space. We, however, showed that in general it seems to be an excellent 
approximation to replace the expectation values of these operators squared by the product 
of expectation values of the individual operators. Only for the very special case of f2 = 2 we 
found that some caution has to prevail, though a perturbative expansion has been already 
proposed (Eq. (|5^)) to improve this approximation when needed. 

Concerning the numerical results we found that SCRPA yields for this model case excel- 
lent results (besides Q = 2, see above). For instance we found that groundstate as well as 
excited energies are always close but consistently above the exact values. We also calculated 
the occupation numbers from the proposed form of the single particle propagator and found 
that they are closest to the exact values in comparison with other proposed approximate 
forms for the occupation numbers. Somewhat as a surprise comes the fact that the so called 
number operator method yields results not better than the quasiboson approximation. We 
give reasons which may back that this is in fact a generic feature. One should say, however, 
that the numerical differences for the occupation numbers using the different methods are, 
at least for the model considered, not very pronounced. 

We also should mention that it is not very difficult to obtain good results for the Lipkin 
model in incorporating groundstate correlations in one way or the other. However, at com- 
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parable numerical complexity, the SCRPA equations do at least equally well, if not better, 
than any other theory on the market. In this respect we refer the reader to our earlier 
study of ref. 0]. A more severe test would be to apply the present SCRPA scheme to other 
more realistic models like for example the multilevel pairing model for which, in the super- 
fluid phase, the number operator approximation is not anymore valid. Such studies shall be 
presented in future work. 
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Table Captions 

Table 1: Correlation energy Ec^yr^. as a function of the interaction strength x- 
Table 2: Excitation energy £■ as a function of the interaction strength x- 
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Table 1 



X 


RPA 


RRPA 


SCRPA 


exact 


0.00 


0.00000 


0.00000 


0.00000 


0.00000 


0.05 


-0.00072 


-0.00067 


-0.00067 


-0.00067 


0.10 


-0.00289 


-0.00270 


-0.00270 


-0.00270 


0.15 


-0.00653 


-0.00608 


-0.00608 


-0.00608 


0.20 


-0.01167 


-0.01083 


-0.01083 


-0.01084 


0.25 


-0.01836 


-0.01698 


-0.01698 


-0.01700 


0.30 


-0.02666 


-0.02455 


-0.02456 


-0.02458 


0.35 


-0.03665 


-0.03358 


-0.03359 


-0.03364 


0.40 


-0.04846 


-0.04409 


-0.04411 


-0.04420 


0.45 


-0.06221 


-0.05612 


-0.05618 


-0.05634 


0.50 


-0.07809 


-0.06972 


-0.06985 


-0.07011 


0.55 


-0.09635 


-0.08490 


-0.08517 


-0.08560 


0.60 


-0.11731 


-0.10166 


-0.10220 


-0.10288 


0.65 


-0.14142 


-0.11996 


-0.12101 


-0.12206 


0.70 


-0.16932 


-0.13963 


-0.14165 


-0.14325 


0.75 


-0.20199 


-0.16029 


-0.16418 


-0.16660 


0.80 


-0.24103 


-0.18117 


-0.18864 


-0.19224 


0.85 


-0.28936 


-0.20064 


-0.21506 


-0.22035 


0.90 


-0.35353 


-0.21552 


-0.24344 


-0.25111 


0.95 


-0.45504 


-0.21994 


-0.27375 


-0.28474 


1.00 




-0.20449 


-0.30596 


-0.32145 


1.05 




-0.15752 


-0.34001 


-0.36151 


1.10 




-0.07081 


-0.37584 


-0.40517 


1.15 




0.05370 


-0.41340 


-0.45271 


1.20 




0.20349 


-0.45263 


-0.50440 


1.25 




0.36207 


-0.49349 


-0.56051 


1.30 




0.51542 


-0.53594 


-0.62129 
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Table 2 



X 


RPA 


RRPA 


SCRPA 


exact 


0.00 


1.00000 


1.00000 


1.00000 


1.00000 


0.05 


0.99875 


0.99875 


0.99894 


0.99894 


0.10 


0.99499 


0.99499 


0.99577 


0.99577 


0.15 


0.98869 


0.98870 


0.99048 


0.99048 


0.20 


0.97980 


0.97986 


0.98308 


0.98308 


0.25 


0.96825 


0.96840 


0.97359 


0.97356 


0.30 


0.95394 


0.95426 


0.96202 


0.96194 


0.35 


0.93675 


0.93737 


0.94840 


0.94821 


0.40 


0.91652 


0.91762 


0.93290 


0.93240 


0.45 


0.89303 


0.89491 


0.91539 


0.91450 


0.50 


0.86603 


0.86908 


0.89605 


0.89455 


0.55 


0.83516 


0.83998 


0.87500 


0.87258 


0.60 


0.80000 


0.80744 


0.85242 


0.84862 


0.65 


0.75993 


0.77126 


0.82852 


0.82275 


0.70 


0.71414 


0.73126 


0.80358 


0.79503 


0.75 


0.66144 


0.68728 


0.77795 


0.76555 


0.80 


0.60000 


0.63925 


0.75200 


0.73444 


0.85 


0.52678 


0.58729 


0.72616 


0.70184 


0.90 


0.43589 


0.53194 


0.70088 


0.66793 


0.95 


0.31225 


0.47437 


0.67658 


0.63290 


1.00 




0.41663 


0.65362 


0.59701 


1.05 




0.36150 


0.63226 


0.56050 


1.10 




0.31175 


0.61269 


0.52369 


1.15 




0.26907 


0.59499 


0.48690 


1.20 




0.23372 


0.57914 


0.45046 


1.25 




0.20494 


0.56507 


0.41472 


1.30 




0.18157 


0.55267 


0.38001 
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Figure Captions 

Figure 1: Selfconsistent mean field potential for quantum fluctuations. 
Figure 2: The mass operator of the single particle Dyson equation represented in terms 
of the full two body T-matrix. 

Figure 3: Feynman graph representing E^^^. 

Figure 4: Groundstate graph for the mass operator G^^. 

Figre 5: The ratio r = —\J {Jq) / (Jo) for the fixed interaction strength x — 
— l)/£ = lasa function of Q,. 

Figure 6: Correlation energy E^orr vs. the interaction strength x, using the exact solu- 
tions (full squares), the RPA (dashed hne), the RRPA (small dots) and the SCRPA (full 
hne) . 

Figure 7: Excitation energy E vs. the interaction strength with the same convention 
of Fig. 6. 

Figure 8: Occupation numbers Q/2 + (Jq) as a function of 
Figure 9: Occupation numbers Q/2 + ( Jq) as a function of Q. for fixed x = 1- 
Figure 10: RPA components of the wave function Y/ X as function of the interaction 
strength. The hues follow the same convention as in Fig. 6. 
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